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The spin of conduction electrons decays due to the combined effect of spin-orbit couphng 
and momentum scattering. The spin-orbit couphng couples the spin to the electron momen- 
tum, that is randomized by momentum scattering off of impurities and phonons. Seen from 
the perspective of the electron spin, the spin-orbit coupling gives a spin precession, while 
momentum scattering makes this precession randomly fluctuating, both in magnitude and 
orientation. 

The speciflc mechanisms for the spin relaxation of conduction electrons were proposed by 
Elliott [Ij and Yafet [2], for conductors with a center of inversion symmetry, and by D'yakonov 
and Perel' [3J, for conductors without an inversion center. In p-doped semiconductors there 
is in play another spin relaxation mechanism, due to Bir, Aronov, and Pikus |,4J. As this has 
a rather limited validity we do not describe it here. More details can be found in reviews 



Before we discuss the two main mechanisms, we introduce a toy model that captures 
the relevant physics of spin relaxation without resorting explicitly to quantum mechanics: 
the electron spin in a randomly fluctuating magnetic field. We will flnd certain universal 
qualitative features of the spin relaxation and dephasing in physically important situations. 

The next part of this review covers the experimental as well as computational status of 
the fleld, discussing the spin relaxation in semiconductors under varying conditions such as 
temperature and doping density. 



0.1 Toy model: the electron spin in a fluctuating magnetic fleld. 

Consider an electron spin S (or the corresponding magnetic moment) in the presence of 
an external time-independent magnetic fleld Bq = Bqz giving rise to the Larmor precession 
frequency ujq = Uqz, and a fluctuating time- dependent fleld B(t) giving the Larmor frequency 
see Fig. [1] We assume that the fleld fluctuates about zero and is correlated on the 
time scale of t^. 



Here a and /3 denote the cartesian coordinates and the overline denotes averaging over 
different random realizations B(t). We will see later that such fluctuating flelds arise quite 
naturally in the context of the electron spins in solids. 

The following description applies equally to the classical magnetic moment described by 
the vector S as well as to the quantum mechanical spin whose expectation value is S. Writing 

out the torque equation, S = a; x S, we get the following equations of motion: 



[SlillTllHlE]. 





(1) 



Sy 



-UoSy + UJy{t)Sz - UJz{t)S, 

uqSx - uJx{t)S^ + uj^{t)Sx, 

Ux{t)Sy - UJy{t)Sx- 



(2) 
(3) 
(4) 



These equations are valid for one speciflc realization of oj{t). Our goal is to flnd instead 

effective equations for the time evolution of the average spin, S(t), given the ensemble of 
Larmor frequencies ijj{t). 



0. 1 . TOY MODEL: THE ELECTRON SPIN IN A FL UCTUATING MAGNETIC FIELD. 3 




Figure 1: Electron spin precesses about the static Bq field along z. The randomly fluctuating 
magnetic fleld B(t) causes spin relaxation and spin dephasing. 

It is convenient to introduce the complex "rotating" spins 5*^ and Larmor frequencies u± 
in the (x, y) plane: 

S+ = Sx -\- iSy, = Sx iSy, (5) 
0J+ = cox + icoy, uj^ = ujx — i(jJy (6) 

The inverse relations are 

S. = + 5, = ^(5+-5_), (7) 

ujx = + . • 

The equations of motion for the spin set (5"+, S^, S^) are, 

5'+ = iuoS^ + iujzS^ — iu+Sz, (9) 
S- = —iuoS- — iUzS- + iu^Sz, (10) 
Sz = -{l/2i){uj+S--uj-S+). (11) 

In the absence of the fluctuating flelds the spin rotates in the complex plane anticlockwise 
(for uq > 0), while clockwise. 

The precession about Bq can be factored out by applying the ansatz, Q: 

S± = s±(t)e^'"«*. (12) 

Indeed, it is straightforward to flnd the time evolution of the set (s+, s_, = S^): 

s+ = iUzS+ — iu-^-Sze~^^°\ (13) 
s_ = -iuzs^ + iuj-Sze''^°\ (14) 
Sz = -(l/2i) (w+s_e-^'^«* -cj_s+e^'^o*) . (15) 

The penalty for transforming into this "rotating frame" is the appearance of the phase factors 



^This is analogous to going to the interaction picture when dealing with a quantum mechanical problem of that type. 
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The solutions of Eqs. [131 EH and [151 can be written in terms of the integral equations, 

s+{t) = s+{0)+i [ dt'u,{t')s+{t')~i [ dt'uj+{t')s,{t')e-''^'''\ (16) 
Jo Jo 

s_{t) = s_{0)-i [ dt'uj,{t')s_{t') + i f dt'co4t')s,it')e''^°'' , (17) 
Jo Jo 
1 /■* r 

Sz{t) = s,(0)--/ dt' u;+(t')s-(t')e-*""*' -w-(^')s+(0e'""*' ■ fl^ 



2« _„ 

We should now substitute the above solutions back into Eqs. [131 [HI and [El The corre- 
sponding expressions become rather lengthy, so we demonstrate the procedure on the s+ 
component only. We get. 



s+{t) = iLo,{t)s+{0) - Lo,{t) / dt'Lo,{t')s+{t') +uj,{t) / dt'Lo+{t')s,{t')e 



-iu+{t)e-''^°'s,{0) + ^e-^^''*a;+(t) / dt' \u+{t')s-{t')e-'''°'' - U-{t')s+{t')e''^fih 

The reader is encouraged to write the analogous equations for (that for s_ is easy to write 
since s_ = s^). 

We now make two approximations. First, we assume that the fluctuating field is rather 
weak and stay in the second order in co^ This allows us to factorize the averaging over the 
statistical realizations of the field, 



u{t)u{t')s{t') ^ u{t)u{t') s{t') (20) 

as the spin changes only weakly over the time scale, r^, of the changes of the fluctuating 
fields. This approximation is called the Born approximation, alluding to the analogy with 
the second-order time-dependent perturbation theory in quantum mechanics. Going beyond 
the Born approximation one would need to execute complicated averaging schemes of the 
product in Eq. [20l since s(t) in general depends on u(t' < t). 

As the second assumption, we consider a "coarse-grained" time evolution, meaning that we 
are effectively averaging s{t) over the time scale of the correlation time Tc] we are interested 
in times t much greater than t^. That allows us to approximate. 



dt' u{t)u{t') s{t') ^ dt'u{t)u{t') s{t), (21) 



since the correlation function u{t)u{t') is significant in the time interval of \t — t'\ ^ Tc only. 
The above approximation makes clear that the spin s(t) is the representative coarse-grained 
(running-averaged) spin of the time interval (t — Tc,t). Equation [21] is a realization of the 
Markov approximation. The physical meaning is that the spin s varies only slowly on the 
time scale of Tc over which the correlation of the fluctuating fields is significant. We then 
need to restrict ourselves to the time scales t larger than the correlation time r^. In effect, 
we will see that in this approximation the rate of change of the spin at a given time depends 
on the spin at that time, not on the previous history of the spin. 

^More precisely, we assume that |a;(t)|Tc <C 1, so that the spin does not fully process about the fluctuating field before the 
field makes a random change. 
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Applying the Born-Markov approximation to Eq. [19] we obtain for the average spin s+ 
the following time evolution equation^ 



Jo Jo 



-iu;+(t)e-*'^o*s^(0) + -e 



-iojot 



dt' 



Using the rules of Eqs. [T] the above simplifies to 



j cit'e-(*-*')/^=s+(t) + -e-'"°*y dt' 



Since we consider the times t ^ Tc, we can approximate 



(23) 



Similarly, 



1 =F ic^oTc 
1 + uir^ 



(24) 



(25) 



The imaginary parts induce the precession of s±, which is equivalent to shifting (renormal- 
izing) the Larmor frequency Uq. The relative change of the frequency is (cuTc)^ which is 
assumed much smaller than one by our Born approximation. We thus keep the real parts 
only and obtain, 



1 



2 1 + U^T^ 



(26) 



Using the same procedure (or simply using s_ = s^) we would arrive for the analogous 
equation for s_: 



— — ^ 1 



Similarly, 



(27) 



(28) 



For t/ie rest of the section we omit the overline on the symbols for the spins, so that S 
will mean the average spin. Returning back to our rest frame of the spins rotating with 
frequency cuq, we get 



1 r, 



2 1 + WqV^ 
1 



2 1 + u;oV2 



1 + a;oV2 



(29) 

(30) 
(31) 



^The initial values of the spin, s(0), arc fixed and not affected by averaging. 
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Finally, going back to Sx and Sy-. 

Sx = -oJoSy - uj^TcSx - ^ o ^lSx (32) 

Sy = a;o^x - a;gTe5y - ^ ^ (^1)'^;/, (33) 

i -|- UJqT^ 

= -(^ + ^)— ^5.. (34) 

We can give the above equation a more conventional form, by introducing two types of the 
spin decay times. First, we define the spin relaxation time Ti by, 

and the spin dephasing times T2 by 

¥- = ^^^ + r^j^' ^^^^ 

W = + '''' 



We then write: 



S 

Sx = -i^oSy - — , (38) 



2a; 



Sy = a;oS.-|^, (39) 



(40) 



Our fiuctuating field is effectively at infinite temperature, at which the average value for the 
spin in a magnetic field is zero. A more general spin dynamics is 

Sx = -uoSy-^, (41) 

-I-2X 

Sy = ^oSx-^, (42) 



where Sqz is the equilibrium value of the spin in the presence of the static magnetic field of 
the Larmor frequency at the temperature at which the environmental fields giving rise 
to u}{t) are in equilibrium. The above equations are called the Block equations. 

The spin components Sx and Sy, which are perpendicular to the applieed static field Bq, 
decay exponentially on the time scales of T2X and T2y, respectively. These times are termed 
spin dephasing times, as they describe the loss of the phase of the spin components perpen- 
dicular to the static field Bq. They are also often called transverse times, for that reason. 
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The time Ti is termed the spin relaxation time, as it describes the (thermal) relaxation of 
the spin to the equilibrium. During the spin relaxation in a static magnetic field the energy 
is exchanged with the environment. In the language of statistical physics, the relaxation 
process establishes the Boltzmann probability distribution for the system. Similarly, de- 
phasing estabhshes the "random phases" postulate that says that there is no correlation 
(coherence) among the degenerate states, such as the two transverse spin orientations; in 
thermal equilibrium such states add incoherently. 

For the sake of discussion consider an isotropic system in which 

^^^ = ^,=^. (44) 

If the static magnetic field is weak, (jJqTc <^ 1, the three times are equal: 

Ti.=T2, = r2, = ^. (45) 

There is no difference between the spin relaxation and spin dephasing. It is at first sight 
surprising that the spin relaxation time is inversely proportional to the correlation time. The 
more random the external field appears, the less the spin decays. We will explain this fact 
below by the phenomenon of motional narrowing. 

In the opposite limit of large Larmor frequency, ujoTc ^ 1, the spin relaxation rate van- 
ishes, 

1 ^ 1 



y - , ^0, (46) 

il Wq Tc 

while the spin dephasing time is given by what is called secular broadening, 

^ ~ ^^c- (47) 

If secular broadening is absent, the leading term in the dephasing time will be, as in the 
relaxation. 

In this limit the spin dephasing rate is proportional to the correlation rate, not to the 
correlation time. 

In the cases in which there is no clear distinction between Ti and T2, we often use the 
symbol 

r, = Ti = T2, (49) 

to describe the generic spin relaxation. 



Motional narrowing 



The surprising fact that, at low magnetic fields, the spin relaxation rate is proportional 

to the correlation time of the fluctuating fleld (as opposed to its inverse), is explained by 
motional narrowing. Consider the spin transverse to an applied magnetic field and assume 
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Figure 2: Electron spin precesses about the static Bq field along z perpendicular to this 
page. The spatially fiuctuating magnetic field Biz causes reversible spin spin dephasing. 

that the field has a single magnitude, but can randomly switch directions, between up to 
down, leading to a random precession of the spin clock and anticlockwise. In effect, the spin 
phase executes a random walk. A single step takes the time t^, the correlation time of the 
fluctuating fleld. After n steps, that is, after the time t = nr, the standard deviation of the 
phase will be 6(j) = {ujTc)y/n, the well known result for a random walk. We call the spin 
dephasing time the time it takes for 5(f) ^ 1. This happens after the time Ts = T^j {utc)'^, or 
Ts = l/{u'^Tc), which is the result we obtained earlier from the Born-Markov approximation, 
Eq. m 



0.1.1 Reversible dephasing, spin ensemble, random walk in inhomogeneous 
fields. 

Our previous calculation was carried out for a single spin in a fluctuating magnetic fleld. Af- 
ter the decay of the spin components, the information about the original spin is irreversibly 
lost as we have no information on the actual history of the fluctuating fleld. Such an irre- 
versible loss of spin is often termed spin decoherence. We will see that spin dephasing can 
be reversible. We will also see that a simple exponential decay, of the type exp(— t/rs), is 
not a rule. 



Reversible spin dephasing: spin ensemble in spatially random magnetic field 

There are physically relevant cases in which the decay of spin is reversible. A typical example 
is an ensemble of localized spins, each precessing about a static local magnetic fleld Bq + 
Bi, with varying static components Bi (of zero average) giving rise to random precession 
frequencies uji. See Fig. [2l Another important example is that of the conduction electrons 
in noncentrosymmetric crystals, such as GaAs, in which the spin-orbit coupling acts as a 
momentum dependent magnetic field. The spins of the electrons of different momenta precess 
with different frequencies. We are interested in the total spin as the sum of the individual 
spins. 



Consider the external field along z direction, and take the fiuctuating frequencies from 
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the Gaussian distribution: 

Piu,) = ^=i=e^-^/^'-^ (50) 

V iTTOW^ 

with zero mean and Scul variance. Denote the in-plane spin of the electron a by and 
Sy. This spin precesses with the frequency UQ + uf, leading to the time evolution for the 
rotating spins 

Sl{t) = S^{t) ± tS^yit) = ^^(0)e^^'^°*e±*'^f*. (51) 

Suppose at t = all the spins are lined up, that is, 5*^(0) = S'±(0). The total spin S±(t) is 
the sum, 

/oo 
dcoiP{uJi)e^'^^\ (52) 

Evaluating the Gaussian integral we get 

S±{t) = 5±(0)e±'"o*e-^"?*'/l (53) 

The in-plane component vanishes after the time of about l/6ui, but this dephasing of the 
spin is reversible, since each individual spin preserves the memory of the initial state. The 
disappearance of the spin is purely due to the statistical averaging over an ensemble in which 
the individual spins have, after certain time, random phases. This spin decay is not a simple 
exponential, but rather Gaussian. 



Spin echo 



We have seen that there are irreversible and reversible effects both present in spin dephasing. 
It turns out that the reversible effects can be separated out by the phenomenon of the spin 
echo. Figure [3] explains this mechanism in detail. Suppose all the localized spins in our 
ensemble point in one direction at time t = 0. At a later time, t = T^, the spins dephase due 
to the inhomogeneities of the precession frequencies and the total spin is small. We apply a 
short pulse of an external magnetic field, the so called tt pulse, that rotates the spins along 
the axis parallel to the original spin direction, mapping the spins as (5*2;, Sy) — )■ {—Sx, Sy). At 
that moment the spins will still be dephased, but the one that is the fastest is now the last, 
and the one that is the slowest, appears as the first. At the time t = 2T.,^ all the spins catch 
on, producing a large spin signal along the original spin direction. In reality this signal will 
be weaker than that at t = due to the presence of irreversible processes, by exp(— 2T^/T2). 
Important, reversible processes are not counted in T2. 



Spin random walk in inhomogeneous magnetic field 



Another interesting situation appears when we consider the possibility that the spin diffuses 
through a region of an inhomogeneous precession frequency^] We can model this situation 
on the system of spins in one dimension, each spin jumping in a time r left or right. Suppose 
the precession frequencies vary in the x direction as 

ijj{x) = Uq + u'x. (54) 



*This discussions applies to nuclear spins as well. 

^The inhomogeneity could be due to the magnetic field or to the g-factor. 
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Figure 3: Initially all the spins point up. Due to random precession frequencies the spins 
soon point in different directions and the average spin vanishes. Applying a vr-pulse rotating 
the spins along the vertical axis makes the fastest spin the slowest, and vice versa. After 
the fastest ones catch up again with the slowest, the original value of the average spin is 
restored. Any reduction from the original value is due to irreversible processes, signal. 

Here u' is the gradient of u. As the electron diffuses, it precesses with the precession 
frequency u{t) = uj[x(t)], given by the position of the electron x{t) at time t. The time 
evolution for an individual spin is 



We are interested in the averaged quantities over many realizations of the random walk. We 
also consider the time scales much greater than the individual random walks steps r. 

The accumulated phase after N steps or t = Nt time is the sum 



Here 5j = ±1 is a random variable representing the random step left or right. The variance 
of (pN is 



According to the central limit theorem, (pjsf is distributed normally with the above variance: 




(55) 



0JV = xi + a;2 + ... +xn = N6i + {N - + ... + Sn- 



(56) 




N-l 



(57) 




(58) 



We can now transform the dynamical equation into the ensemble averaging: 




(59) 



where we denoted the diffusivity hj D = t/2, for the unit length step. 
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Figure 4: The electron performs a random walk. Its spin precesses by the inhomogeneous 
magnetic field along the x axis. The average spin dephases to zero with time. 



0.1.2 Quantum mechanical description 



The Born-Markov approximation to the dynamics of a system coupled to an external envi- 
ronment can be cast in the quantum mechanical language. We refer the reader to Ref. f8] 
for the derivation. The physics of this derivation is the same as what we did in Sec. 10. 11 
Only the formalism is different. 



We consider the system described by the Hamiltonian 

H{t) = Ho + V{t), 



(60) 



in which Hq is our system per se and V{t) is the time-dependent random fluctuating fleld, 
of zero average and correlation time Tc. 



V{t) = V^t) = 0, V{t)Vit') ~ e-l*-*'l/"= 
The system is fully described by the density matrix p. 



(61) 



The transformation to the rotating frame is equivalent to going to the interaction picture 
in quantum mechanics: 



Piit) 
Vj{t) 



iHot/h 



V(t)e 



-iHot/h 



(62) 
(63) 



Performing the operations as outlined in Sec. 10.11 for the classical model, we arrive at the 
effective time evolution for the density of state operator, 



dpijt) 
dt 



dt' 



Vj{t), [Vi{t'),pi{t)] 



(64) 



The above equation is called the Master equation and is the starting equation in many 
important problems in which a quantum system is in contact with a reservoir. 
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Figure 5: D'yakonov-Perel' mechanism. The electron starts with the spin up. As it moves, its 
spin precesses about the axis corresponding to the electron velocity. Phonons and impurities 
change the velocity, making the spin to precess about a different axis (and with different 
speed). During the scattering event the spin is preserved. 

0.1.3 Spin relaxation of conduction electrons 

We will consider nonmagnetic conductors in zero or weak magnetic fields so that the spin 
dephasing and spin relaxation times are equal, T2 = Ti = Tg. The formula for the spin 
relaxation time, 

1/ts = uj^Tc, (65) 

that was derived above for the spin in a fluctuating magnetic field, applies in a semiquan- 
titative sense (that is, it gives an order of magnitude estimates and useful trends) to the 
conduction electron spins as well. We analyze below the D'yakonov-Perel' |3] and the EUiott- 
Yafet[Tll2] mechanisms. 

The D'yakonov-Perel' mechanism is at play in solids lacking a center of spatial inversion 
symmetry. The most prominent example is the semiconductor GaAs. In such solids the 
spin-orbit coupling is manifested as some effective magnetic field — the spin- orbit field — that 
depends on the electron momentum. Electrons in different momentum states feel different 
spin-orbit fields, so that the spin precesses with a given Larmor frequency, until the electron 
is scattered into another momentum state. See Fig. |5l As the electron momentum changes 
on the time scale r of the momentum relaxation time, the net effect of the momentum 
scattering on the spin is to produce random fluctuations of the Larmor frequencies. We have 
motional narrowing. Since these frequencies are correlated by Tc = r, we arrive at 

— = ^sV, (66) 

for the spin relaxation time. The magnitude of Ugo is the measure of the strength of the spin- 
orbit coupling. The spin relaxation rate is directly proportional to the momentum relaxation 
time: the more the electron scatters, the less its spin dephases. 

The Elliott- Yafet mechanism works in systems with and without a center of inversion. 
It relies on spin-flip momentum scattering. The spin-flip amplitudes are due to spin-orbit 
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Figure 6: Elliott-Yafet mechanism. The electron starts with the spin up. As it scatters off 
of impurities and phonons, the spin can also flip due to spin-orbit coupling. Between the 
scattering the spin is preserved. After, say, a million scattering events, the spin will be down. 

coupling, while the momentum scattering is due to the presence of impurities (that also 
contribute to the spin-orbit coupling), phonons, rough boundaries, or whatever is capable 
of randomizing the electron momentum. During the scattering events, the spin is preserved. 
See Fig. O How do we account for such a scenario with our toy model? Consider an electron 
scattering off of an impurity with a spin flip. This spin flip can be viewed as a precession 
that occurs during the time of the interaction of the electron with the impurity. Let us say 
that the scattering takes the time of Xp/vp, where vp is the electron velocity and Xp is the 
electron wave-length at the Fermi level (considering that it is greater or at most equal to the 
the size a of the impurity — otherwise we could equally put a/vp). Then the precession angle 
{p = Uso{Xp/vp), with Uso denoting the spin-orbit coupling induced precession frequency. 
Let us compare this with the angle of precession, ut, in the motional narrowing model, 
Eq. [651 We see that the spin-flip can be described by the effective precession frequency 
^ = <^so{Xp/vpT). We then obtain for the spin relaxation rate 



Here Ego = f^so and e = {hkp)vp/2 is the Fermi energy. For the the Elliott-Yafet mechanism 
holds that the more the electron scatters, the more the spin dephases. 

Below we discuss the two mechanisms in more detail, providing the formalisms for their 
investigation. 



0.2 The D'yakonov-Perel' mechanism 

D'yakonov and Perel' [3] considered solids without a center of inversion symmetry, such as 
GaAs or InAs. The mechanism also works for electrons at surfaces and interfaces. In such 
solids the presence of spin-orbit coupling induces the spin-orbit fields which give rise to the 
spin precession. Momentum relaxation then causes the random walk of the spin phases. 




(67) 
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Spin-orbit field 

In solids without a center of inversion symmetry, spin-orbit coupling splits the electron 
energies: 

£k,t 7^ (68) 
Only the Kramers degeneracy is left, due to time reversal symmetry: 

ek,t = e-Ki- (69) 

This energy splitting at a given momentum k is conveniently described by the spin-orbit field 
Q, giving a Zeeman-like (but momentum dependent) energy contribution to the electronic 
states, described by the additional Hamiltonian (to the usual band structure): 

H, = ■ a. (70) 

The time reversal symmetry requires that the spin-orbit field is an odd function of the 
momentum: 

f^k = -f^-k. (71) 

The representative example of a spin-orbit field is the Bychkov-Rashba [10] (a^ij) and Dres- 
selhaus pT^TS] (7/5) couplings in 2d electron gases formed at the zinc-blende heterostructures 
grown along [001] [8j: 

Hi = {aBR + lD)<^xky - {aBR - lD)<^ykx- (72) 

The axes are x = [110] and y = [110]. The spin-orbit field 

hQk = 2 [{aBR + lD)ky, -{aBR - IdiK] ■ (73) 

This field has the C2v symmetry, reflecting the structural symmetry of the zinc-blend inter- 
faces (such as GaAs/GaAlAs) with the principal axes along [110] and [110]. 



Kinetic equation for the spin 

Let Sk be the electron spin in the momentum state k. The time evolution of the spin is then 
described by 

^ - flk X Sk = - V Wkk' (sk - si,) . (74) 

The left hand side gives the full time derivative dsi^/dt, which is due to the explicit change 
of the spin, its direction, and the change of the spin due to the change of the momentum k. 
The right-hand side is the change of the spin at k due to the spin-preserving scattering from 
and to that state. The scattering rate between the two momentum states k and k' is Wkk'- 

There are two time scales in the problem. First, momentum scattering, which occurs on 
the time scale of the momentum relaxation time r, makes the spins in different momentum 
states equal. Second, this uniform spin decays on the time scale of the spin relaxation time 
Tg which we need to find. We assume that r ^ r^; this assumption is well satisfied in real 
systems. In principle we could go directly to our model of the electron spin in a fluctuating 
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magnetic field, with the role of the random Larmor precession playing by 17, identifying 
Tc = T. However, it is instructive to see how this two time-scale problem is solved directly. 

We separate the fast and slow components of the spins, 

Sk = s + (Ik) = 0. (75) 

The symbol (...) denotes averaging over different momenta . Our goal is to find the effective 
equation for the time evolution of s, which is the actual spin averaged over k. The fast 
component, ^i^, decays on the time scale of r to the value given by the instantaneous value 
of s. Our goal is to find the effective equation for the time evolution of s, which is the 
actual spin averaged over k. The fast component, ^j^., decays on the time scale of r to the 
quasistatic value given by the instantaneous value of s. 

Upon substituting Eq. [75] into the kinetic equation Eq. [71] and averaging over k, we 
obtain the equation of motion for the averaged spin 

s=(f7kX^k,) (76) 

using that (fik) = 0. We need to find ^j^.. Since s is hardly changing on the time scales 
relevant to ^j^., we write, 

- X s - flk X = - E ^kk' (Ck - Ik') • (77) 

k' 

We make the following assumption: 

< 1. (78) 

That is, we assume that the precession is slow on the time scale of the momentum relaxation 
time (see our note on the Born approximation in the toy model section Sec. 10.11 For 
simplicity, we make the relaxation time approximation to model the time evolution of the 
fast component ^k- 

^k - f7k X s - flk X Ik = - — • (79) 

T 

From the condition of the quasistatic behavior, d^^^/dt = 0, we get up to the first order in 
Qt the following solution for the quasistatic |k- 

Ik = r(f7kxs). (80) 

The above is a realization of coarse graining, in which we effectively average the spin evolution 
over the time scales of r. 

Substituting the quasistatic value of |k from Eq. [80] to the time evolution equation for s, 
Eq. [76] we find 

s = r(rik X (rik X s)). (81) 
Using the vector product identities we finally obtain for the individual spin components a, 

Sa = {^ka^ki3)rsp - {fllr)sa. (82) 

These equations describe the effective time evolution of the electron spin in the presence of 
the momentum dependent Larmor precession ilk- 
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For the specific model of tlie zinc-blend heterostructure with the C2V spin-orbitc field, Eq. 
[73| we obtain the spin dephasing dynamics: 

where 

1 [aBR ^ipf 1 1 2 
— = — 2 — ' ~ = ~' y^^i 

and 

- = -|^£fc(«L + 7l))- (85) 

The up (down) sign is for the Sx (sy). The spin relaxation is anisotropic. The maximum 
anisotropy is for the case of equal magnitudes of the Bychkov-Rashba and Dresselhaus in- 
teractions, ubr = ±7d- In this case one of the spin components does not decay. The 
component of the spin relaxes roughly twice faster than the in-plane components. 



0.2.1 The persistent spin helix 



Let us consider the case of asR = Id = A/2. According to Eq. IHHthe spin component Sx 
does not decay, while the decay rates of Sy and are the same: 

Sy = -2sy/Ts, s^ = -2sJts. (86) 

It turns out that a particular nonuniform superposition of Sy and can exhibit no decay as 
well. This superposition has been termed persistent spin helix [T3] . 



Let us assume that the spin is no longer uniform, so that the kinetic equation contains 
the spin gradient as well, due to the quasiclassical change of the electronic positions: 

^ - f^k X Sk + ^ ■ Vk = - V W^kk' (Sk - s;,) . (87) 
ot or ^-^ 

k' 

A running spin wave, 

Sk(r) = Skc''^"; Sk = Sk(q), (88) 

then evolves according to, 

^ - rik X Sk + z (q ■ Vk) Sk = - ^ W^kk' (sk - Sk) . (89) 

k' 

We again separate the fast and slow spins, 

Sk = s + Ik, (Ik) = 0. (90) 
For the dynamics of the slow part we get, 

Sk=(f^kX|k)-^((q-Vk)lk)- (91) 



^The decay of that component would be due to higher-order (such as cubic) terms in the spin-orbit fields. 



0.2. THE D 'YAKONOV-PEREL ' MECHANISM 



17 



Proceeding as in the previous section, assuming that s is stationary on the time scales 
relevant to ^, we can write, 

- f^k X s - flk X + « (q ■ Vk) s + i (q ■ Vk) = - 5^ W^kk' (Ik - Ik') • (92) 

k' 

We now work with the following assumptions: 

fir < 1, g£ < 1, (93) 

where i = gr is the mean free path. We thus assume that the precession is slow on the time 
scale of the momentum relaxation time, as well as (this is new here) the electronic motion 
is diffusive on the scale of the wavelength of the spin wave. In the momentum relaxation 
approximation, also considering the leading terms according to the conditions Eq. [931 "we 

get 

ik-f^kXS + z(q-Vk)s = -l (94) 

r 

In the steady-state, corresponding to a given s(t), the solution is 

Ik = T (f^k X s) - zr (q ■ Vk) s. (95) 
Substituting to the time evolution equation for s, Eq. [93l we find 

s = -2zr((q ■ Vk) (flk x s)) + r{n^ x {il^ x s)) - r((q ■ Vk)')s. (96) 
Using the vector product identities and introducing the diffusivity 

D = (vL)r, (97) 
where a denote the cartesian coordinates (we assume an isotropic system), we finally obtain 

Sa = -2iTea0yqs{Vk5^kp)s^ + {flka^k0)rSp - {fliT)sa - Dq^Sa. (98) 

For our specific case of 

fik = (Afcj„0,0), X = aBR + l3D. (99) 
we find that Sx decays only via diffusion: 

Sx = -Dqh,. (100) 
The spin dephasing is ineffective. This is an expected result. 

More interesting behavior is found for the two remaining spin components, Sy and s^. 
These two spin components, transverse to the spin-orbit field, are coupled: 

Sy = +2i^qyDs, - (rfi2 + Dq^) Sy, (101) 
Trt\ 

= -2i—qyDsy - (rfi2 + Dq^) s,, (102) 

and we denoted = {fll). Interestingly, this set of coupled differential equations has only 
decaying solutions. This is best seen by looking at the rotating spins: 

s+ = Sy + isz, = Sy — isz, (103) 
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Figure 7: The persistent spin helix is a wave of circularly polarized spin. The sense of 
polarization, clock or counterclockwise, depends on the relative sign of the Dresselhaus and 
the Bychkov-Rashba spin-orbit coupling. 



whose time evolutions are uncoupled. 



-(n'T + Dq' - '^qyD ] (104) 
2mA 



- ( fiV + ^ -^qyD ] s_ (105) 



Considering that 



fiV = X^kl)r = (^) {v^y)r = (^) D, (106) 
we find that the decay of s+ vanishes for the wave vector 

poTT mX m , ^ . / _x 

gf"" = = ^ («^« + /^«) • (107) 

The abbreviation PSH stands for the persistent spin helix, which is the spin wave described 
by s+ at this particular wave vector. See Fig. [71 While individually both Sy and decay 
at a generic wave vector (and also in the uniform case, g = 0), the spin helix they form 
does not decay in the approximation of the linear spin-orbit field. The spin wave rotating in 
the opposite sense, s_, on the other hand, decays. And vice versa for q^^^ = —mX/h. The 
persistent spin helix was observed in the spin grating experiment 



0.3 The Elliott- Yafet mechanism 



Elliott fTj was first to recognize the role of the intrinsic spin-orbit coupling — that coming 
from the host ions — on spin relaxation. Yafet [2] significanly extended this theory to properly 
treat electron-phonon spin-fiip scattering. The Elliott- Yafet mechanism dominates the spin 
relaxation of conduction electrons in elemental metals and semiconductors, the systems 
with space inversion symmetry. In systems lacking this symmetry, the mechanism appears 
competes with the D'yakonov-Perel' one; the dominance of one over the other depends on 
the material in question and specific conditions, such as temperature and doping. 
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Suppose a nonequilibrium spin accumulates in a nonmagnetic degenerate conductor. The 
spin accumulation is given as the difference between the chemical potentials for the spin up 
and the spin down electrons. Let us denote the corresponding potentials by /J,^ and //j,. The 
nonequilibrium electron occupation function for the spin A is 



1 



g/3(£k-MA) + 1 



dSy- 



Here, 



(108) 



(109) 



describes the equilibrium degenerate electronic system of the Fermi energy ep- The electron 
density for the spin A is 



riA = y" d£gs{e)fl{£) + j ^^^(^) 



df{e) 
de 



(110) 



Here Qsi^s) is the electronic density of states, defined per unit volume and per spin, at the 
fermi level: 



9s 



E 



The total electron density n is 



n 



I deg{e)f{e). 



(Ill) 



(112) 



We assume that the spin accumulation does not charge the system, that is, the charge 
neutrality is preserved n^ + n^ = n. This condition is well satisfied in metals and degenerate 
semiconductors that we consider. We then get 



The spin density is. 



where is the spin quasichemical potential, — (A*t ~ 



(113) 
(114) 



The spin relaxation time Ti is defined by the decay law. 



ds 
di 



duf dn 



dt dt ■■n-^nt--7^--ys7^- (us) 

Here W^j, is the net number of transitions per unit time from the spin 4- to t- Similarly, W_i_-^ 
expresses the rate of spin flips from f to J,. In the degenerate conductors the spin decay is 
directly proportional to the decay of the spin accumulation /x^: 

dus _ 
'dt ~ ¥i 



(116) 



0.3.1 The electron-impurity scattering 



We need to distinguish the cases of the impurity or host induced spin-orbit couphng. Al- 
tahugh the formulas for the calculation of Ti look similar in the two cases, they are never- 
theless conceptually different. 
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Spin-orbit coupling by the impurity 



If the spin-orbit couphng comes from the impurity potential (in this case the coupling is 
often termed extrisnic) the spin-flip scattering is due to that potential. 

The number of transitions per unit time from the spin up to the spin down states is 

^n = ^Yl ^k'n't.kni - W^k<,k'n't- (117) 
kn k'n' 

The rate is given by the spin-flip events from down to up minus the ones from up to down. 
We use the Fermi golden rule to write out the individual scattering rates: 



tiknj. 



/kn (1 — /k'n) l^knt.k'naP"^ (^^k'n' — ^kn) , 



knt,k'n'^ 



2tt 

-^fk'n' (1 — /kn) |f4;'n'4.,kntl ^ (^k'n' — £kn) ■ 



118) 

;ii9) 



Substituting for the occupation numbers the l inearized expression, Eq. I108[ and using the 
definition of the spin relaxation of Eq. 1115^ we obtain for the spin relaxation rate the 
expression, 



^-y-Z.I^'^-t-k'n' 



4-1 



kk' 



S (^k'n' — ^kn) • 



If we define the spin relaxation for the individual momentum state k as 



7^ - ^ |t^knt,k'n'4 

-ilk n '—f 



^ (f^k'n' — £^kn) 



(120) 



(121) 



which has a straightforward interpretation as the spin-flip rate by the elastic impurity scat- 
tering to all the possible states k', we get 



1 



1 



:i22) 



as the average of the individual scattering rates over the electronic Fermi surface. 



Spin-orbit coupling by the host lattice 



If the spin-orbit coupling comes from the host lattice (in this case the coupling is often 
termed intrinsic) the spin-flip scattering is due to the admixture of the Pauli spin up and 
spin down states in the Bloch eigenstates. 

How do the Bloch states actually look like in the presence of spin-orbit coupling? Elliott 
showed that the Bloch states corresponding to a generic lattice wave vector k and band n 
can be written as pLj, 

*k,nt(r) = [akn(r)| t) + &kn(r)U)]e*'^^ (123) 
^k,n^(r) = [al,„(r)|i)-r,„(r)U)]e^'^-. (124) 
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The states | t) and | 4-) are the usual Pauh spinors. We can select the two states such that 
l^knl ~ 1 while \bkn\ <C 1, due to the weak spin orbit coupling; this justifies calling the two 
above states "spin up" (ft) and "spin down" (JJ,)- In fact, to "prepare" the states for the 
calculation of the spin relaxation, they need to satisfy, 

{k,n\\a,\k,n\') = \5xy, (125) 

with A =f|', Jj.. That is, the two states should diagonalize the spin matrix Sz (or whatever 
spin direction one is interested in). 

The Bloch states of Eqs. 11231 and 11241 allow for a spin fiip even if the impurity does not 
induce a spin-orbit coupling. Indeed, the matrix element 

{k,nii\U\k,nij.) r^ab, (126) 

is in general non zero due to the spin admixture. The spin fiip probability is proportional to 
the spin admixture probability. This quantity is crucial in estimating the spin relaxation 
in the Elliott- Yafet mechanism. The spin relaxation time Ti in this case can be calculated 
using the formula Eq. 11201 with 

f4;nt,k'na = f4;nfr,k'n'4]., (127) 

given by Eq. 11261 A useful rule of thumb for estimating the spin relaxation time in this case 
is 

1 ^ i^, (128) 

where the averaging of the spin admixture probabilities 6^ is performed over the Fermi surface 
(or the relevant energy scales of the problem); Tp is the spin-conserving momentum relaxation 
time. We stress that b is obtained from the states prepared according to Eq. 11251 



0.3.2 The electron-phonon scattering 



The spin-fiip due to the scattering of the electrons off of phonons involves the intrinsic spin- 
orbit potential. The elect ron Bloch states are the ones given by the Eqs. 11231 and 11241 
prepared according to Eq. 11251 

The net number of transitions per unit time from the spin up to the spin down states is 

W^ti ~ W^knt,qJ';k'n'4. + W^knt;k'n'4,,q!/ — Wk'n'4,,qi/;knt ~ Wk'n'ijknt.q;'- (129) 

kn k'n' qy 

We introduced the rates of the spin flip transitions accompanied by the phonon absorption 
and emissions as follows. The net transition rate from the single electron state {k'n' \.) to 
the electron state |k?7, '\) while the phonon of momentum q and polarization i/ is emitted, is 

W^knt,q!^;k'n'4. = | ^knt,qi^;k'n'4. | /k'n'J, (1 — /knf) ^ (^kn " ^k'n' + f^qu) ■ (130) 

Similarly, the net transition rate from the single electron state \k'n' \.) to the electron state 
\kn '\) while the phonon of momentum q and polarization u is absorbed, is 

2vr 2 

WkntM'n'l,qu = | A^knt;k'n'4.qi/ 1 fk'n'l (1 — /knt) ^ (^kn — ^k'n' ~ ^qi/) • (131) 
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The same way are defined the remaining two rates, VFk'n'4„qi/;knt spin flip from kn 

to k'n' I with the phonon emission, and Wk'n'4.;kntqi') for the phonon absorption: 

W^k'n'4.,qz/;knt = | ^k'n'4.,qi/;knt I /knf (1 ~ fk'n'i) ^ {^k'n' — ^kn + ^<^(iu) , (132) 

W^k'n'i;kntq;^ = |^k'n'4„q;/;knt,qi^l /knf ~ /k'n'J,) (s^k'n' " ^^kn " ^qi/) • (133) 

The calculation of the spin relaxation due to the electron-phonon scattering is rather 
involved and we cite here only the flnal result for degenerate conductors: 

^= T-Y.Y.Y.(-—) Jm W.-.:.-{knt\VV\k'u'li)\' (134) 

-^1 ^ 5's f;-^ ^ V ^kn / 2A/M Wk'-k,!. 

kn k'n' u ^ ' ' 

X { [n-qu - fu'n' + l] (^kn " £^k'n' " ^qi/) + ["-qi. + /k'n'] ^ (^kn " £^k'n' + ^q(^|5) 

The electronic bands are described by the energies £kn of the state with momentum k and 
band index n. Phonon frequencies are Wqjy, for the phonon of momentum q and polarization 
z/; similarly for the phonon polarization vector eq^. We further denoted by M the atomic 
mass, by the number of atoms in the lattice, and by Vl^ the gradient of the electron- 
lattice ion potential. The equilibrium phonon occupation numbers are denoted as nqj,, given 
as ^ 

n^u = nM = (136) 

The electronic states |kncr) are normalized to the whole space. 

Two types of processes contribute to the phonon induced spin flips. First, what we call 
the Elli ott p rocesses, are the Elliott-type of spin flips in which the Bloch states given by Eq. 
I123l and 11241 scatter by the scalar part of the gradient of the electron-ion potential V . Second, 
what we call the Yafet processes, are the spin flips due to the gradient of the spin-orbit part 
of the electron-ion potential. These two processes are typically of similar order of magnitude 
and have to be added coherently in order to obtain Ti. 

Figure [H] shows the experimental data of the spin relaxation in intrinsic (nondegenerate) 
silicon, obtained by the spin resonance |8] and the spin injection [17J experiments. The 
calculation based on the Elliott- Yafet mechanism of the phonon induced spin flips reproduces 
the experiments [T5] . 



The Yafet relation 

The expected temperature dependence of the phonon-induced spin flips in degenerate con- 
ductors, is, according to the Elliott- Yafet mechanism, 1/Ti ~ T at high temperatures (start- 
ing roughly at a fraction of the Debye temperature Tjj) and 1/Ti ~ at low temper- 
atures, in analogy with the conventional spin-conserving electron-phonon scattering. The 
high-temperature dependence originates from the linear increase of the phonon occupation 
numbers n with increasing temperature: nqj, ~ T, as ksT ^ ujf^y. At T > To are all the 
phonons excited. The low temperature dependence of the spin- conserving scattering follows 
from setting the relevant phonon energy scale to ksT. The matrix element 

(kn^ |V\/|kV^) (137) 
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Figure 8: Spin relaxation in silicon. Phonon- induced spin relaxation in silicon results in an 
approximate power law [15]. Shown are the experimental data from spin resonance [16] 
and spin injection [17] experiments, and a calculation based on the Elliott- Yafet mechanism 
115]. 



which gives the 1/Ti ~ dependence. However, Yafet showed [2] that for the spin- flip 
matrix element the space inversion symmetry modifies the momentum dependence to 

(kn ^ |Vr|kV ^) ~ (138) 

so that 

1 ~ T^ (139) 

J 1 

instead of the expected T^. Since the same temperature dependence holds for the phonon- 
induced electrical resistance p(T), we can write 

^-piT), (140) 
J 1 

known as the Yafet relation. 

The relation Eq. 11381 holds if both the Elliott and Yafet processes are taken into account. 
Individually, they would lead to a linear dependence on q. The quantum mechanical inter- 
ference between these two processes thus significantly reduces the spin-flip probability at low 
momenta q. An example is shown in Fig. |9l The Elliott and Yafet processes would indi- 
vidually give much stronger spin relaxation than is observed. Their destructive interference 
can modify Ti by orders of magnitude. 



0.4 Results based on kinetic spin Bloch equation approach 



It was shown by Wu et at. from a full microscopic kinetic-spin-Bloch-equation approach [9] 
that the single-particle approach is inadequate in accounting for the spin relaxation/dephasing 
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Figure 9: Interference between the Elliott and Yafet processes. Individually the Elliott and 
Yafet phonon- induced spin relaxation processes give spin relaxation orders of magnitude 
stronger that the total one. This example is from the calculation of Ti in silicon |15] . 

both in the time |18l [191 l20l [21] [22] and in the space [231 l2ll [25] domains. The momentum 
dependence of the effective magnetic field (the D'yakonov-Perel' term) and the momentum 
dependence of the spin diffusion rate along the spacial gradient |23| or even the random spin- 
orbit coupling [26j all serve as inhomogeneous broadenings [T9l I2U]. It was pointed out that 
in the presence of inhomogeneous broadening, any scattering, including the carrier-carrier 
Coulomb scattering, can cause an irreversible spin relaxation/dephasing [19]. Moreover, be- 
sides the spin relaxation/dephasing channel the scattering provides, it also gives rise to the 
counter effect to the inhomogeneous broadening. The scattering tends to drive carriers to a 
more homogeneous state and therefore suppresses the inhomogeneous broadening. Finally, 
this approach is valid in both strong and weak scattering limits and also can be used to study 
systems far away from the equilibrium, thanks to the inclusion of the Coulomb scattering. 

In the following, we present the main results based on kinetic-spin-Bloch-equation ap- 
proach. We first briefly introduce the kinetic spin Bloch equations. Then we review the 
results of the spin relaxation/dephasing in the time and space domains respectively. 



0.4.1 Kinetic spin Bloch equations 

By using the nonequilibrium Green function method with gradient expression as well as the 
generalized Kadanoff-Baym Ansatz [27], we construct the kinetic spin Bloch equations as 
follows: 



Here p]f_{r,t) are the density matrices of electrons with momentum k at position r and 
time t. The off-diagonal elements of pk represent the correlations between the conduction 
and valence bands, different subbands (in confined structures) and different spin states. 
Pk(r,t)|^j- are the driving terms from the external electric field. The coherent terms in Eq. 



Pk(r,t) 



Pk(r,t)|d^ + Pk(r,t)| 



dif 




(141) 
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(1141 p pklcoh are composed of the energy spectrum, magnetic field and effective magnetic 
field from the D'yakonov-Perel' term, and the Coulomb Hartree-Fock terms. The diffusion 
terms Pk(r,^)ldif come from the spacial gradient. The scattering terms Pk(r,i)|scat include 
the spin-flip and spin conserving electron-electron, electron-phonon and electron-impurity 
scatterings. The spin-flip terms correspond to the Elliot- Yafet and/or Bir-Aronov-Pikus 
mechanisms. Detailed expressions of these terms in the kinetic spin Bloch equations depend 
on the band structures, doping situations and dimensionalities [9J and can be found in the lit- 
erature for different cases, such as intrinsic quantum wells [18j, n-type quantum wells without 
[211 EHl EHl pi Ip^ and with [22l El [33l IS] electric field, p-type quantum wells [Ml ESI ET] , 
quantum wires [581 139] , quantum dots [^ and bulk materia ls [^Lj in the spacial uniform case 
and quantum wells in spacial non- uniform case [23l [2H |42l [43] . By numerically solving the 
kinetic spin Bloch equations with all the scattering explicitly included, one is able to obtain 
the time evolution and/or spacial distribution of the density matrices, and hence all the mea- 
surable quantities, such as mobility, diffusion constant, optical relaxation/dephasing time, 
spin relaxation/dephasing time, spin diffusion length, as well as hot-electron temperature, 
can be determined from the theory without any fitting parameters. 



0.4.2 Spin relaxation/dephasing 

In this subsection we present the main understandings of the spin relaxation/dephasing added 
to the literature from the kinetic-spin-Bloch-equation approach. We focus on three related 
issues: (i) The importance of the Coulomb interaction to the spin relaxation/dephasing; (ii) 
Spin dynamics far away from the equilibrium; and (iii) Qualitatively different behaviors from 
those wildly adopted in the literature. 

First we address the effect of the Coulomb interaction. Based on the single-particle 
approach, it has been long believed that the Coulomb scattering is irrelevant to the spin 
relaxation/dephasing [H]. It was first pointed out by Wu and Ning [19] that in the presence 
of inhomogeneous broadening in spin precession, i.e., the spin precession frequencies are It- 
dependent, any scattering, including the spin-conserving scattering, can cause irreversible 
spin dephasing. This inhomogeneous broadening can come from the energy-dependent g- 
factor fT9], the D'yakonov-Perel' term 20j, the random spin-orbit coupling [26], and even 
the momentum dependence of the spin diffusion rate along the spacial gradient [23]. Wu and 
Ning first showed that with the energy-dependent (^-factor as an inhomogeneous broadening, 
the Coulomb scattering can lead to irreversible spin dephasing [19]. In [001] -grown n-doped 
quantum wells, the importance of the Coulomb scattering for spin relaxation/dephasing was 
jroved by Glazov and Ivchenko [15] by using perturbation theory and by Weng and Wu 
21] from the kinetic-spin-Bloch-equation approach. In a temperature-dependent study of 
;he spin dephasing in [001]-oriented n-doped quantum wells, Leyland et al. experimentally 
verified the effects of the electron-electron Coulomb scattering by closely measuring the 
momentum scattering rate from the mobility [16]. By showing the momentum relaxation 
rate obtained from the mobility cannot give the correct spin relaxation rate, they showed 
the difference comes from the Coulomb scattering. Later Zhou et al. even predicted a peak 
from the Coulomb scattering in the temperature dependence of the spin relaxation time in 
a high-mobility low-density n-doped (001) quantum well [29j. This was later demonstrated 
by Ruan et al. experimentally [17] . 

Figure [10] shows the temperature dependence of the spin relaxation time of a 7.5 nm 
GaAs/Alo.4Gao.6As quantum well at different electron and impurity densities [29]. For this 
small well width, only the lowest subband is needed in the calculation. It is shown in the 



26 




50 100 150 200 250 300 
T(K) 



Figure 10: Spin relaxation time r vs. the temperature T with well width a = 7.5 nm 
and electron density n being (a) 4 x 10^*^ cm~^, (b) 1 x 10^^ cm~^, and (c) 2 x 10^^ cm~^, 
respectively. Solid curves with triangles: impurity density = ra; solid curves with dots: 
Hi = O.ln; solid curves with circles: rii = 0; dashed curves with dots: Ui = O.ln and no 
Coulomb scattering. From Zhou et al. [29]. 



figure that when the electron-impurity scattering is dominant, the spin relaxation time de- 
creases with increasing temperature monotonically. This is in good agreement with the 
experimental findings [48] and a nice agreement of the theory and the experimental data 
from 20 to 300 K is given in Ref. [29]. However, it is shown that for sample with high mobil- 
ity, i.e., low impurity density, when the electron density is low enough, there is a peak at low 
temperature. This peak, located around the Fermi temperature of electrons T|. = Ep/kB, 
is identified to be solely due to the Coulomb scattering [221 SS]- It disappears when the 
Coulomb scattering is switched off, as shown by the dashed curves in the figure. This peak 
also disappears at high impurity densities. It is also noted in Fig. [TOT c) that for electrons of 
high density so that Tp is high enough and the contribution from the electron-longitudinal 
optical-phonon scattering becomes marked, the peak disappears even for sample with no im- 
purity and the spin relaxation time increases with temperature monotonically. The physics 
leading to the peak is due to the crossover of the Coulomb scattering from the degener- 
ate to the non-degenerate limit. At T < Tp, electrons are in the degenerate limit and the 
electron-electron scattering rate l/rge oc T^. At T > Tp, 1/ree oc T~^ [121 ED]- Therefore, at 
low electron density so that Tp is low enough and the electron-acoustic phonon scattering is 
very weak comparing with the electron-electron Coulomb scattering, the Coulomb scattering 
is dominant for high mobility sample. Hence the different temperature dependence of the 
Coulomb scattering leads to the peak. It is noted that the peak is just a feature of the 
crossover from the degenerate to the non-degenerate limit. The location of the peak also 
depends on the strength of the inhomogeneous broadening. When the inhomogeneous broad- 
ening depends on momentum linearly, the peak tends to appear at the Fermi temperature. 
A similar peak was predicted in the electron spin relaxation in p-type GaAs quantum well 
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and the hole spin relaxation in (001) strained asymmetric Si/SiGe quantum well, where the 
electron and hole spin relaxation times both show a peak at the hole Fermi temperature Tp 
[371 ES]. When the inhomogeneous broadening depends on momentum cubically, the peak 
tends to shift to a lower temperature. It was predicted that a peak in the temperature de- 
pendence of the electron spin relaxation time appears at a temperature in the range of (T|./4, 
Tp/2) in the intrinsic bulk GaAs [H] and a peak in the temperature dependence of the hole 
spin relaxation time at Tp/2 in p-type Ge/SiGe quantum well [37j. Ruan et al. demonstrated 
the peak experimentally in a high-mobility low-density GaAs/Alo.ssGao.esAs heterostructure 
and showed a peak appears at Tp/2 in the spin relaxation time versus temperature curve 

For larger well width, the situation may become different in the non-degenerate limit. 
Weng and Wu calculated the spin relaxation/dephasing for (001) GaAs quantum wells with 
larger well width and high mobility, by including the multi-subband effect [28]. It is shown 
that for small/large well width so that the linear /cubic Dresselhaus term is dominant, the 
spin relaxation/dephasing time increases/decreases with the temperature. This is because 
with the increase of temperature, both the inhomogeneous broadening and the scattering 
get enhanced. The relative importance of these two competing effects is different when the 
linear/cubic term is dominant [28j. Jiang and Wu further introduced strain to change the 
relative importance of the linear and cubic D'yakonov-Perel' terms and showed the different 
temperature dependences of the spin relaxation time [51]. This prediction has been real- 
ized experimentally by Holleitner et al. where they showed that in n-type two-dimensional 
InGaAs channels, when the linear D'yakonov-Perel' term is suppressed, the spin relaxation 
time decreases with temperature monotonically [5_2j. Another interesting prediction related 
to the multi-subband effect is related to the effect of the inter-subband Coulomb scattering. 
From the calculation Weng and Wu found out that although the inhomogeneous broaden- 
ing from the higher subband of the (001) quantum well is much larger, due to the strong 
inter-subband Coulomb scattering, the spin relaxation times of the lowest two subbands are 
identical [28]. This prediction has later been verified experimentally by Zhang et al., who 
studied the spin dynamics in a single-barrier heterostructure by time-resolved Kerr rotation 
[53] . By applying a gate voltage, they effectively manipulated the confinement of the second 
subband and the measured spin relaxation times of the first and second subbands are almost 
identical at large gate voltage. Lii et al. showed that due to the Coulomb scattering, T2 = Tg* 
in (001) GaAs quantum wells for a wide temperature and density regime [51]. It was also 
pointed out by Lii et al. that in the strong (weak) scattering limit, introducing the Coulomb 
scattering will always lead to a faster (slower) spin relaxation/dephasing [35] . 

Another important effect from the Coulomb interaction to the spin relaxation/dephasing 
comes from the Coulomb Hartree-Fock contribution in the coherent terms of the kinetic spin 
Bloch equations. Weng and Wu [21] first pointed out that at a high spin polarization, the 
Hartree-Fock term serves as an effective magnetic field along the z axis which blocks the 
spin precession. As a result, the spin relaxation/dephasing time increases dramatically with 
the spin polarization. They further pointed out that the spin relaxation/dephasing time 
decreases with temperature at high spin polarization in quantum well with small well width, 
which is in contrast to the situation with small spin polarizations. These predictions have 
been verified experimentally by Stich et al. in an n-type (001) GaAs quantum well with high 
mobility [551 ESj. By changing the intensity of the circularly polarized lasers, Stich et al. 
measured the spin dephasing time in a high mobility n-type GaAs quantum well as a function 
of initial spin polarization. Indeed they observed an increase of the spin dephasing time with 
the increased spin polarization, and the theoretical calculation based on the kinetic spin 
Bloch equations nicely reproduced the experimental findings when the Hartree-Fock term 
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Figure 11: (a) The spin dephasing times as a function of initial spin polarization for constant, 
low excitation density and variable polarization degree of the pump beam. The measured 
spin dephasing times are compared to calculations with and without the Hartree-Fock (HF) 
term, showing its importance, (b) The spin dephasing times measured and calculated for 
constant, high excitation density and variable polarization degree. From Stich et al. [56] . 



was included [55] . It was also shown that when the Hartree-Fock term is removed, one does 
not see any increase of the spin dephasing time. Later, they further improved the experiment 
by replacing the circular-polarized laser pumping with the elliptic polarized laser pumping. 
By doing so, they were able to vary the spin polarization without changing the carrier density. 
Figure. [TT] shows the measured spin dephasing times as function of initial spin polarization 
under two fixed pumping intensities, together with the theoretical calculations with and 
without the Coulomb Hartree-Fock term. Again the spin dephasing time increases with 
the initial spin polarization as predicted and the theoretical calculations with the Hartree- 
Fock term are in good agreement with the experimental data [56]. Moreover, Stich et al. 
also confirmed the prediction of the temperature dependences of the spin dephasing time 
at low and high spin polarizations [56]. Figure ITW a) shows the measured temperature 
dependences of the spin dephasing time at different initial spin polarizations. As predicted, 
the spin dephasing time increases with increasing temperature at small spin polarization but 
decreases at large spin polarization. The theoretical calculations also nicely reproduced the 
experimental data. The hot-electron temperatures in the calculation were taken from the 
experiment [Fig. [T2t^b)]. The effective magnetic field from the Hartree-Fock term has been 
measured by Zhang et al. from the sign switch of the Kerr signal and the phase reversal 
of Larmor precessions with a bias voltage in a GaAs heterostructure [57]. Korn et al. [58] 
also estimated the average effect by applying an external magnetic field in the Faraday 
configuration, as shown in Fig. fTST a) for the same sample reported above [551 ES]. They 
compared the spin dephasing times of both large and small spin polarizations as function of 
external magnetic field. Due to the effective magnetic field from the Hartree-Fock term, the 
spin relaxation times are different under small external magnetic field but become identical 
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Figure 12: (a) Spin dephasing time as a function of sample temperature, for different ini- 
tial spin polarizations. The measured data points are represented by solid points, while 
the calculated data are represented by lines of the same colour, (b) Electron temperature 
determined from intensity-dependent photoilluminance measurements as a function of the 
nominal sample temperature, for different pump beam fluence and initial spin polarization, 
under experimental conditions corresponding to the measurements shown in (a). The mea- 
sured data points are represented by solid points, while the curves serve as guide to the eye. 
From Stich et al. 



when the magnetic field becomes large enough. From the merging point, they estimated 
the mean value of the effective magnetic field is below 0.4 T. They further showed that this 
effective magnetic field from the Hartree-Fock term cannot be compensated by the external 
magnetic field, because it does not break the time-reversal symmetry and is therefore not 
a genuine magnetic field, as said above. This can be seen from Fig. [TST b) that the spin 
relaxation time at large spin polarization shows identical external magnetic field dependences 
when the magnetic field is parallel or antiparallel to the growth direction. 
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Figure 13: (a) Spin dephasing times as a function of an external magnetic field perpendicular 
to the quantum well plane for small and large initial spin polarization, (b) Same as (a) for 
large initial spin polarization and both polarities of the external magnetic field. From Korn 
et al.[58j 

We now turn to discuss the spin relaxation/dephasing far away from the equilibrium. In 
fact, the spin relaxation/dephasing of high spin polarization addressed above is one of the 
cases far away from the equilibrium. Another case is the spin dynamics in the presence of a 
high in-plane electric field. The spin dynamics in the presence of a high in-plane electric field 
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was first studied by Weng et al. [22] in GaAs quantum well with only the lowest subband by 
solving the kinetic spin Bloch equations. To avoid the "runaway" effect [59] , the electric field 
was calculated upto 1 kV/cm. Then Weng and Wu further introduced the second subband 
into the model and the in-plane electric field was increased upto 3 kV/cm [2B]- Zhang et 
al. included L valley and the electric field was further increased upto 7 kV/cm [32J. The 
effect of in-plane electric field to the spin relaxation in system with strain was investigated 
by Jiang and Wu [51]. Zhou et al. also investigated the electric-field effect at low lattice 
temperatures [29] . 



The in-plane electric field leads to two effects: i) It shifts the center-of-mass of electrons 
to krf = m*Vrf = m*/iE with /i representing the mobility, which further induces an effective 
magnetic field via the D'yakonov-Perel' term [22]. ii) The in-plane electric field also leads to 
the hot-electron effect [6U]. The first effect induces a spin precession even in the absence of 
any external magnetic field and the spin precession frequency changes with the direction of 
the electric field in the presence of an external magnetic field [221 ESj- The second effect en- 
hances both the inhomogeneous broadening and the scattering, two competing effects leading 
to rich electric- field dependence of the spin relaxation/dephasing and thus spin manipulation 
[22lE8l[32l[5llE9lE3lil]. 

Finally we address some issues of which the kinetic-spin-Bloch-equation approach gives 
qualitatively different predictions from those widely used in the literature. These issues in- 
clude the Bir-Aronov-Pikus mechanism, the Elliot- Yafet mechanism and some density/temperature 
dependences of the spin relaxation/dephasing time. 
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Figure 14: Ratio of the spin relaxation time due to the Bir-Aronov-Pikus mechanism to 
that due to the D'yakonov-Perel' mechanism, tbap/^dp, as function of temperature and hole 
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satisfying tbap/tdp = 1- Note the smaller the ratio rBAp/T"DP is, the more important the 
Bir-Aronov-Pikus mechanism becomes. The yellow solid curves indicate the cases satisfying 
9^^[A^Ljj(i) + N^^{2)]/dfj_i^Nh = 0.1. In the regime above the yellow curve the multi-liole- 
subband effect becomes significant. From Zhou et al. [36] . 



It has long been believed in the literature that for electron relaxation/dephasing, the 
Bir-Aronov-Pikus mechanism is dominant at low temperature in j9-type samples and has im- 
portant contribution to intrinsic sample with high photo-excitation [61, 62, 63, 64, 65, 66. 67j. 
These conclusion was made based on the single-particle Fermi golden rule. Zhou and We re- 
examined the problem using the kinetic-spin-Bloch-equation approach [30] . They pointed out 
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that the Pauh blocking was overlooked in the Fermi Golden rule approach. When electrons 
are in the non-degenerate limit, the results calculated from the Fermi Golden rule approach 
are valid. However, at low temperature, electrons can be degenerate and the Pauli blocking 
becomes very important. As a result, the previous approaches always overestimated the 
importance of the Bir-Aronov-Pikus mechanism at low temperature. Moreover, the previous 
single-particle theories underestimated the contribution of the D'yakonov-Perel' mechanism 
by neglecting the Coulomb scattering. Both made the Bir-Aronov-Pikus mechanism dom- 
inate the spin relaxation/dephasing at low temperature. Later, Zhou et al. performed a 
thorough investigation of electron spin relaxation in p-type (001) GaAs quantum wells by 
varying impurity, hole and photo-excited electron densities over a wide range of values [36J, 
under the idea that very high impurity density and very low photo-excited electron den- 
sity may effectively suppress the importance of the D'yakonov-Perel' mechanism and the 
Pauli blocking. Then the relative importance of the Bir-Aronov-Pikus and D'yakonov-Perel' 
mechanisms may be reversed. This indeed happens as shown in the phase-diagram-like pic- 
ture in Fig. [m where the relative importance of the Bir-Aronov-Pikus and D'yakonov-Perel' 
mechanisms is plotted as function of hole density and temperature at low and high impurity 
densities and photo-excitation densities. For the situation of high hole density they even 
included multi-hole subbands as well as the light hole band. It is interesting to see from 
the figures that at relatively high photo-excitations, the Bir-Aronov-Pikus mechanism be- 
comes more important than the D'yakonov-Perel' mechanism only at high hole densities and 
high temperatures (around hole Fermi temperature) when the impurity is very low [zero in 
Fig. [Mf a)]. Impurities can suppress the D'yakonov-Perel' mechanism and hence enhance 
the relative importance of the Bir-Aronov-Pikus mechanism. As a result, the temperature 
regime is extended, ranging from the hole Fermi temperature to the electron Fermi temper- 
ature for high hole density. When the photo-excitation is weak so that the Pauli blocking is 
less important, the temperature regime where the Bir-Aronov-Pikus mechanism is important 
becomes wider compared to the high excitation case. In particular, if the impurity density is 
high enough and the photo-excitation is so low that the electron Fermi temperature is below 
the lowest temperature of the investigation, the Bir-Aronov-Pikus mechanism can dominate 
the whole temperature regime of the investigation at sufficiently high hole density, as shown 
in Fig. [T^ d). The corresponding spin relaxation times of each mechanism under high or low 
impurity and photo-excitation densities are demonstrated in Fig. [TH They also discussed the 
density dependences of spin relaxation with some intriguing properties related to the high 
hole subbands [36j. The predicted Pauli-blo eking effect in the Bir-Aronov-Pikus mechanism 
has been partially demonstrated experimentally by Yang et al. [68] They showed by increas- 
ing the pumping density, the temperature dependence of the spin dephasing time deviates 
from the one from the Bir-Aronov-Pikus mechanism and the peaks at high excitations agree 
well with those predicted by Zhou and Wu |30j . 

Another widely accepted but incorrect conclusion is related to the Elliot- Yafet mecha- 
nism. It is widely accepted in the literature that the Elliot- Yafet mechanism dominates spin 
relaxation in n-type bulk III-V semiconductor at low temperature, while the D'yakonov- 
Perel' mechanism is important at high temperature |69l El [71 El EO] . Jiang and Wu pointed 
out that the previous understanding are based on the formula that can only be used in 
the nondegenerate limit. Moreover, the momentum relaxation rates are calculated via the 
approximated formula for mobility [69]. By performing an accurate calculation via the 
kinetic-spin-Bloch-equation approach, they showed that the Elliot- Yafet mechanism is not 
important in III-V semiconductors, including even the narrow-band InAs and InSb [3T]. 
Therefore, the D'yakonov-Perel' mechanism is the only spin relaxation mechanism for n- 
type III-V semiconductors in metallic regime. 

Jiang and We have further predicted a peak in the density dependence of the spin relax- 
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Figure 15: Electron spin relaxation times from the calculation via the kinetic-spin-Bloch- 
equation approach (red solid curve) and from the experiment [72\ (blue •) as function of the 
doping density. The green dashed curve shows the results without hot-electron effect. The 
hot-electron temperature used in the computation is plotted as the blue dotted curve (Note 
the scale is on the right hand side of the frame). The chain curve is the calculated spin 
relaxation time with a fixed hot-electron temperature 80 K, and N^.^ = 6 x 10^^ cm~^. From 
Shen [73]. 



ation/dephasing time in n-type III-V semiconductors where the spin relaxation/dephasing 
is limited by the D'yakonov-Perel' mechanism [41] . Previously, the nonmonotonic density 
dependence of spin lifetime was observed in low-temperature (T <5 K) experiments, where 
the localized electrons play a crucial role and the electron system is in the insulating regime 
or around the metal- insulator transition point |7T]. Jiang and Wu found, for the first time, 
that the spin lifetime in metallic regime is also nonmonotonic. Moreover, they pointed out 
that it is a universal behavior for all bulk III-V semiconductors at all temperature where the 
peak is located at Tp ~ T with Tp being the electron Fermi temperature. The underlying 
physics for the nonmonotonic density dependence in metallic regime can be understood as 
following: In the nondegenerate regime, as the distribution is the Boltzmann one, the density 
dependence of the inhomogeneous broadening is marginal. However, the scattering increases 
with the density. Consequently the spin relaxation/dephasing time increases with the den- 
sity. However, in the degenerate regime, due to the Fermi distribution, the inhomogeneous 
broadening increases with the density much faster than the scattering does. As a result, the 
spin relaxation/dephasing time decreases with the density. Similar behavior was also found 
in two-dimensional system j3T| [37], where the underlying physics is similar. The predicted 
peak was later observed by KrauB et al. [T^ as shown in Fig. [T5] where theoretical calculation 
based on the kinetic spin Bloch equations nicely reproduced the experimental data by Shen 



0.4.3 Spin diffusion/transport 

By solving the kinetic spin Bloch equations together with the Poisson equation self-consistently, 
one is able to obtain all the transport properties such as the mobility, charge diffusion length 
and spin diffusion/injection length without any fitting parameter. It was first pointed out by 
Weng and Wu [23] that the drift- diffusion equation approach is inadequate in accounting for 
the spin diffusion/ transport. It is important to include the off-diagonal term between oppo- 
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site spin bands pkfj, studying the spin diffusion/transport. With this term, electron spin 
precesses along the diffusion and therefore k ■ VrPk(i",^) in the diffusion terms pk(i",^)|dif 
provides an additional inhomogeneous broadening. With this additional inhomogeneous 
broadening, any scattering, including the Coulomb scattering, can cause an irreversible spin 
relaxation/dephasing [23j. Unlike the spin precession in the time domain where the inhomo- 
geneous broadening is determined by the effective magnetic field from the D'yakonov-Perel' 
term, h(k), in spin diffusion and transport it is determined by 

n^ = \gfiBB + h{k)\/K, (142) 

provided the diffusion is along the x-axis |12]. Here the magnetic field is in the Voigt configu- 
ration. Therefore, even in the absence of the D'yakonov-Perel' term h(k), the magnetic field 
alone can provide an inhomogeneous broadening and leads to the spin relaxation/dephasing 
in spin diffusion and transport. This was first pointed out by Weng and Wu back to 2002 
[23] and has been realized experimentally by Appelbaum et al. in bulk silicon [711[75], where 
there is no D'yakonov-Perel' spin-orbit coupling due to the center inversion symmetry. Zhang 
and Wu further investigated the spin diffusion and transport in symmetric Si/SiGe quantum 
wells Hg. 

When B = but the D'yakonov-Perel' term is present, then the inhomogeneous broaden- 
ing for spin diffusion and transport is determined by $7^ = h(k)/fca.. In (001) GaAs quantum 
well where the D'yakonov-Perel' term is determined by the Dresselhaus term [76j, the average 
of fik reads (fik) = C{{ky) — (A;^),0,0) with C being a constant. For electrons in quantum 
well, this value is not zero. Therefore, the spacial spin oscillation due to the Dresselhaus 
effective magnetic field survives even at high temperature when the scattering is strong. This 
effect was first predicted by Weng and Wu by showing a spin pulse can oscillate along the 
diffusion in the absence of the magnetic field at very high temperature [23] . Detailed studies 
were carried out later on this effect |25l [771 112] • The spin oscillation without any applied 
magnetic field in the transient spin transport was later observed experimentally by Crooker 
and Smith in strained bulk system [78]. differing from the two-dimensional case, in bulk the 
average of fik from the Dresselhaus term is zero, since (fik) = C{{ky) — (fc^),0,0) = due 
to the symmetry in the y- and 2;-directions. This is consistent with the experimental result 
that there is no spin oscillation for the system without stress. However, when the stress 
is applied, an additional spin-orbit coupling, namely the coupling of electron spins to the 
strain tensor, appears, which is linear in momentum This additional spin-orbit coupling 
also acts as an effective magnetic field. Therefore, once the stress is applied, one can observe 
spacial spin oscillation even when there is no applied magnetic field |78j . 

Cheng and Wu further developed a new numerical scheme to calculate the spin diffu- 
sion/transport in GaAs quantum wells with very high accuracy and speed [12|. It was 
discovered that due to the scattering, especially the Coulomb scattering, T2 = Tg* is valid 
even in the spacial domain. This prediction remains yet to be verified experimentally. More- 
over, as the inhomogeneous broadening in spin diffusion is determined by |h(k)|/A;j. in the 
absence of magnetic field, the period of the spin oscillations along the x-axis is independent 
on the electric field perpendicular to the growth direction of the quantum well [12] , which is 
different from the spin precession rate in the time domain [22]. This is consistent with the 
experimental findings by Beck et al. [79] . 

Cheng et al. applied the kinetic-spin-Bloch-equation approach to study the spin transport 
in the presence of competing Dresselhaus and Rashba fields [SOl] . When the Dresselhaus and 
Bychkov-Rashba [10] terms are both important in semiconductor quantum well, the total 
effective magnetic field can be highly anisotropic and spin dynamics is also highly anisotropic 
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Figure 16: Schematic of the different directions considered for the spin polarizations [(110), 
(110) and (OOl)-axes] and spin diffusion/injection (x-axis). From Cheng et al. [80j. 
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Figure 17: Spin diffusion length Ld (solid curves) and the inverse of the spin oscillation period 
Lq ^ (dashed curves) for identical Dresselhaus and Rashba coupling strengths as functions of 
the injection direction for different spin polarization directions no, z and fii (fii = z x fio, 
i.e., crystal direction [110]) at T = 200 K. It is noted that the scale of the spin oscillation 
period is on the right hand side of the frame. From Cheng et al. [80j. 

in regards to the spin polarization |81] . For some special polarization direction, the spin 
relaxation time is extremely large [8T | '|82 | |83| |8^ . For example, if the coefficients of the 
linear Dresselhaus and Bychkov-Rashba terms are equal to each other in (001) quantum well 
of small well width and the cubic Dresselhaus term is not important, the effective magnetic 
field is along the [110] direction for all electrons. For the spin components perpendicular 
to the [110] direction, this effective magnetic field flips the spin and leads to a finite spin 
relaxation/dephasing time. For spin along the [110] direction, this effective magnetic field can 
not flip it. Therefore, when the spin polarization is along the [110] direction, the Dresselhaus 
and Bychkov-Rashba terms can not cause any spin relaxation/dephasing. When the cubic 
Dresselhaus term is taken into account, the spin dephasing time for spin polarization along 
the [110] direction is finite but still much larger than other directions [S3]. The anisotropy in 
the spin direction is also expected in spin diffusion and transport. When the Dresselhaus and 
Bychkov-Rashba terms are comparable, the spin injection length for the spin polarization 
perpendicular to [110] direction is usually much shorter than that for the spin polarization 
along [110] direction. In the ideal case when there are only the linear Dresselhaus and 
Bychkov-Rashba terms with identical strengths, spin injection length for spin polarization 
parallel to the [110] direction becomes infinity [B21IS3]- This effect has promoted Schliemann 
et al. to propose the nonballistic spin-field-effect transistor |83] . In such a transistor, a gate 
voltage is used to tune the strength of the Bychkov-Rashba term and therefore control the 
spin injection length. However, Cheng et al. pointed out that spin diffusion and transport 
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actually involve both the spin polarization and spin transport directions [80]. The latter 
has long been overlooked in the literature. In the kinetic-spin-Bloch-equation approach, 
this direction corresponds to the spacial gradient in the diffusion term [pk(i", i)|dif] and 
the electric field in the drifting term [pk(r, ^) |dr]- The importance of the spin transport 
direction has not been realized until Cheng et al. pointed out that the spin transport is 
highly anisotropic not only in the sense of the spin polarization direction but also in the 
spin transport direction when the Dresselhaus and Bychkov-Rashba effective magnetic fields 
are comparable [80j. They even predicted that in (001) GaAs quantum well with identical 
linear Dresselhaus and Bychkov-Rashba coupling strengths, the spin injection along [110] 
or [110 [86] can be infinite regardless o/ the direction of the spin polarization. This can 
be easily seen from the inhomogeneous broadening Eq. (11421) which well defines the spin 
diffusion/transport properties. For the spin diffusion/transport in a (001) GaAs quantum 
well with identical Dresselhaus and Bychkov-Rashba strengths (the schematic is shown in 
Fig. [16] with the transport direction chosen along the x-axis), the inhomogeneous broadening 
is given by [80] 



with 9 being the angle between the spin transport direction (x-axis) and [001] crystal direc- 
tion. It can be splitted into two parts: the zeroth-order term (on k) which is always along 
the same direction of no and the second-order term which comes from the cubic Dressel- 
haus term. If the cubic Dresselhaus term is omitted, the effective magnetic fields for all k 
states align along ng (crystal [110]) direction. Therefore, if the spin polarization is along 
rio, there is no spin relaxation even in the presence of scattering since there is no spin pre- 
cession. Nevertheless, it is interesting to see from Eq. fll43p that when 9 = Stt/A, i.e., the 
spin transport is along the [110] direction, fi^ = 2m*/3no is independent on k if the cubic 
Dresselhaus term is neglected. Therefore, in this special spin transport direction, there is no 
inhomogeneous broadening in the spin transport for any spin polarization. The spin injec- 
tion length is therefore infinite regardless of the direction of spin polarization. This result is 
highly counterintuitive, considering that the spin relaxation times for the spin components 
perpendicular to the effective magnetic field are finite in the spacial uniform system. The 
surprisingly contradictory results, i.e., the finite spin relaxation/dephasing time versus the 
infinite spin injection length, are due to the difference in the inhomogeneous broadening in 
spacial uniform and non-uniform systems. For genuine situation, due to the presence of the 
cubic term, the spin injection length is still finite and the maximum spin injection length 
does not happen at the identical Dresselhaus and Bychkov-Rashba coupling strengths, but 
shifted by a small amount due to the cubic term |80J. However, there is strong anisotropy 
in regards to the spin polarization and spin injection direction, as shown in Fig. [T71 This 
predication has not yet been realized experimentally. However, very recent experimental 
findings on spin helix [ 131 [H] have provide strong evidence to support this predication [87j. 

Now we turn to the problem of spin grating. Transient spin grating, whose spin polar- 
ization varies periodically in real space, is excited optically by two non-collinear coherent 
light beams with orthogonal linear polarization [881 [89l [90] [H] . Transient spin grating tech- 
nique can be used to study the spin transport since it can directly probe the decay rate 
of nonuniform spin distributions. Spin diffusion coefficient Dg can be obtained from the 
transient spin grating experiments [HHl [HU [90]. In the literature, the drift- diffusion model 
was employed to extract Dg from the experimental data. With the drift-diffusion model. 
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Figure 18: ! = + r_)/2 and AT = {T+ - r_)/2 vs. q at T = 295 K. Open 
boxes/triangles are the relaxation rates r+/_ calculated from the full kinetic spin Bloch 
equations. Filled/open circles represent F and AF respectively. Noted that the scale for 
AF is on the right hand side of the frame. The solid curves are the fitting to F and AF 
respectively. The dashed curves are guide to eyes. From Weng et al. [9l] . 

the transient spin grating was predicted to decay exponentially with time with a decay 
rate of F^ = Dgcf + l/r^, where q is the wavevector of the spin grating and is the spin 
relaxation time [881 189] . However, this result is not accurate since it neglects the spin pre- 
cession which plays an important role in spin transport as first pointed out by Weng and 
Wu [22]. Indeed, experimental results show that the decay of transient spin grating takes 
a double-exponential form instead of single exponential one [SH [901 [13] • Also the relation 
which relates the spin injection length with the spin diffusion coefficient Ds and the spin 
relaxation time Lg = I^D^Ts from the drift- diffusion model should be checked. In fact, 
if this relaxation is correct, the above prediction of infinite spin injection length at certain 
spin injection direction for any spin polarization in the presence of identical Dresselhaus and 
Bychkov-Rashba coupling strengths cannot be correct. 



Figure 19: Spin relaxation times r± vs. temperature for (a) high-mobility sample with 
q = 0.58 X 10^ cm~^ and (b) low-mobility sample with q = 0.69 x 10^ cm^^. The dots are 
the experiment data from Ref. [90j. From Weng et al. [91 \ . 

Weng et al. studied this problem from the kinetic-spin-Bloch-equation approach [9T]. 
By first solving the kinetic spin Bloch equations analytically by including only the elastic 
scattering, i.e., the electron-impurity scattering, they showed that the transient spin grating 
should decay double exponentially with two decay rates r±. In fact, non of the rates is 
quadratic in q. However, the average of them reads [HI] 




100 150 200 250 100 150 200 250 300 



T(K) 



T (K) 



F = (F+ + F_)/2 = Dg2 + l/r: 



(144) 



with 1/r^ = {1/ts + 1/tsi)/2, which differs from the current widely used formula by replacing 
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the spin decay rate by the average of the in- and out-of-plane relaxation rates. The difference 
of these two decay rates is a hnear function of the wavevector q when q is relatively large: 

AT = cq + d, (145) 

with c and d being two constants. The steady-state spin injection length Lg and spin pre- 
cession period Lq are then [DI] 

Ls = 2Ds/^/\c^-ADs{l/Ti-d)l (146) 
Lo = 2Ds/c. (147) 

They further showed that the above relations Eqs. (11441) and (I145p are valid even including 
the inelastic electron-electron and electron phonon scatterings by solving the full kinetic 
spin Bloch equations, as shown in Fig. [181 A good agreement with the experimental data 
[90] of double exponential decays t± = F^^ are shown in Fig. [191 Finally it was shown 
that the infinite spin injection length predicted by Cheng et al. in the presence of identical 
Dresselhaus and Bychkov-Rashba coupling strengths [HQ] addressed above can exactly be 

obtained from Eq. (11461) as in that special d = and c = I^D^jr, [91]. 

However, ^D^t^ always remains finite unless = oo. Therefore, Eqs. (114411147]) give the 
correct way to extract the spin injection length from the spin grating measurement. 
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